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Abstract 

The singular nature of a non-integer spiral phase plate allows easy manipulation of spatial degrees 
of freedom of photon states. Using two such devices, we have observed very high dimensional 
{D > 3700) spatial entanglement of twin photons generated by spontaneous parametric down- 
conversion. 

PACS numbers: 03.67.Mn, 42.50.Dv 
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FIG. 1: Schematic drawing of a spiral phase plate (SPP) with a step index t = h s {n — no)/X, where 
h s is the step height, n and no are the refractive indices of the SPP and the surrounding medium, 
respectively, and A is the wavelength of the incident light. 

In the early 1990s, it was realized that the quantum state of a photon in a Laguerre- 
Gaussian laser mode LG' is an eigenstate of the orbital angular momentum (OAM) operator, 



with eigenvalue lh (I integer) 
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The mode functions LGf, (/ = — oo . . . oo, p — . . . oo) 



can be used as a basis of an infinite-dimensional Hilbert space, and thus to describe, within 
the paraxial approximation, a photon state represented by an arbitrary spatial wave function. 

According to theory, the twin photons created in spontaneous parametric down-conversion 
(SPDC) are entangled in the full, infinite-dimensional (spatial) Hilbert space; in fact, they 
are maximally entangled Ul^. Several years ago, Mair et al. \^ demonstrated OAM en- 
tanglement in a 3D subspace of that Hilbert space, in an experiment where they used 
integer-OAM analyzers. It has been recently proposed that spatial entanglement can be 
demonstrated in an m/imte-dimensional subspace 6]. This can be achieved by using frac- 
tional-OAM analyzers, in view of the fact that fractional OAM states are coherent superposi- 
tions of an infinite number of LG^ states with integer OAM. Here we report an experimental 
implementation of this idea, employing two half-integer OAM analyzers; we confirm the 
predicted entanglement of fractional OAM and deduce a dimensionality of D > 3700 per 
photon, the entangled space having a dimensionality D 2 . 

The key element in our approach is the use of a spiral phase plate (SPP) 0,0), shown in 
Fig. ^ as part of our analyzer. This plate imprints an azimuth-dependent phase retardation 
on an incident field. The difference in phase retardation A0 at azimuthal angles and 27r is 
equal to A0 = 2ir£, whereby the plate adds an average OAM of £h per photon. By adjusting 
the SPP parameters (see caption of Fig. ^), the SPP index I can be tuned to any desired 
value. If i is non-integer, the output field is not cylindrically symmetric (cf. Fig. Efb)) and 
carries a mixed screw-edge dislocation in the near field, i.e. a singularity with a transverse 
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FIG. 2: Decomposition of the analyzer mode, associated with an SPP with t = 3.5 and a single- 
mode fiber, into LG modes, (a) The contribution per LG mode for mode orders m + n < 7, where 
I = n — m and p = min(n, m). (b) The fidelity of the mode decomposition as a function of number 
of modes used in the decomposition. Each of the data points is calculated for LG modes of mode 
order up to and including the value indicated on the top axis. The solid line is an accurate fit to 
the numerical data, allowing extrapolation. 

component; this dislocation is associated with the radial edge of the SPP Thus, for all 
non-integer £, the SPP is characterized by two parameters, namely the step index £ itself 
and the orientation of the radial edge of the SPP, a. The latter is varied by rotating the 
SPP around the beam axis. 

Our approach involves the use of two fractional-OAM analyzers (one per arm) in an SPDC 
setup. Such an analyzer consists of a non-integer SPP, a single-mode fiber and a detector; 
the orientations of the radial edges of the SPPs are varied (see Fig. Eta) below). Because 
the SPPs have a non-integer step index £, the fractional OAM states provide a natural basis 
t0 d e SCrlb e the Wphoton state Q. Each e.e.eat o f this .action has, can be written 
as an £ and a dependent superposition of an infinite number of LGj, states, extending over 
both I and p. Because this superposition changes when we reorient the SPP, we explore an 
m/zmte-dimensional subspace of the full, (spatial) Hilbert space by systematically varying 
a (at constant step index £). We stress that the contribution of LG^ modes with p ^ to 
this superposition is significant; this is caused by the fact that the dislocation imprinted by 
the SPP on the incident field has a mixed screw-edge character jjj. 

The most convenient implementation of this scheme is by choosing the non-integer part 
C of the step index £ to be half- integer, i.e. £ = Irthe states with a = and a = tc are then 
orthogonal (independent of the integer part of £) [6( . When we place an analyzer containing 
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an SPP with index i in the signal path of a SPDC setup, and one containing an SPP with 
index —I in the idler path, the coincidence probability as a function of the orientations of 
the two SPPs is predicted to be [f| 



viz. a parabolic coincidence fringe with a visibility of unity. The constant C depends on the 
overlap of the radial part of the two-photon state as emanating from the SPDC crystal, with 
the radial part of the fiber mode [6] . This parabolic coincidence fringe, which is a function of 
only the relative orientation of the two SPPs, reveals the entanglement of the twin photons. 

The key to success in the experiment lies in the analyzers or, more specifically, in the two 
SPPs. Our SPPs have been manufactured using a polymer molding technique starting from 
a micro-machined mold, yielding identical plates with I = ±3.48 ± 0.02 at the (degenerate 
SPDC) wavelength of 813 nm. The step of each SPP is not abrupt but in fact a linear ramp 
with an azimuthal width of 6°. Furthermore, the plates show a central 'anomaly' with a 
diameter of about 300 /zm |9(. The rest of the plate has an rms deviation of ~ 15 nm from 
the ideal shape |lfl ]. 

The high-dimensional nature of the i = 3.5 SPP is illustrated in Fig. Efa), which shows 
a modal decomposition for mode order < 7 if the input is a fundamental Gaussian LG[j 
(i.e. the fiber mode). The fidelity of this truncated decomposition is only 43% (Fig. Efb)). 
When truncating at mode order 20 (corresponding to 231 LG modes), the fidelity becomes 
71%. Extrapolation shows that a fidelity of 98% requires about 125000 LG modes. We have 
checked that these results remain essentially unchanged when using the actual topography 
of our SPP plates; in particular, the finite azimuthal width of the step has no significant 
effect. 

The SPDC setup is shown in Fig. |3fa). The pump laser is a 50 mW Krypton-ion laser 
(Spectra Physics 2060) operating at a wavelength of 406.7 nm. Its output is weakly focused 
on a nonlinear crystal (BBO), with beam waist wq = 0.78 ± 0.01 mm. The crystal has a 
thickness L = 1.0 mm and has been cut for nearly collinear type-II SPDC. Before insert- 
ing a fractional-OAM analyzer in each arm, we test the setup by observing polarization 
entanglement; we achieve 95% and 92% fringe visibility in the crystal (H,V) and 45° basis, 
respectively 

We fix the polarizers in the crystal basis and insert a fractional-OAM analyzer with SPP 
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FIG. 3: (a) The SPDC setup for establishing high-dimensional entanglement with two analyzers. 
The SPDC crystal (BBO) is tilted so that the angle between the signal and idler arms is 6°. In 
both arms, one encounters a 0.5 mm-thick wave plate (WP) for walk-off compensation, a telescope 
that images the nonlinear crystal onto the SPP, a Polarcor™ polarizer (P) and a red filter (not 
shown) to eliminate stray pump light. The down-converted photons are focused with microscope 
objectives onto single-mode fibers leading to single-photon-counting modules; coincidence counts 
are measured. Both SPPs can be rotated (a s and on) around the respective beam axes, (b) An 
experimental far-field intensity pattern of the analyzer mode. The white curves indicate the 90% 
iso-intensity lines of the SPDC rings, illustrating that the analyzer mode completely falls in the 
region of twin-photon degeneracy. 

index £ s = +3.48 in the signal path and an analyzer with SPP index £^ = —3.48 in the idler 
path. The SPDC cone-crossings in the far field are found to be much larger than the far- 
field intensity pattern of the analyzer mode, as illustrated in Fig. Of b). The SPDC intensity 
varies about 10% over the SPP diffraction profile. 

The SPPs have to be aligned with great precision. In particular, both the rotation axis 
and the center of each SPP have to be aligned, with a precision of order 1 fim, with the 
core of the corresponding single-mode fiber. Only when this amount of care is invested in 
aligning the optical elements, can one rotate both SPPs to make use of the a s and ctj degrees 
of freedom without losing any signal. 

Our experimental results are shown in Fig. |U four coincidence fringes (obtained for idler 
SPP-settings aj = 0, n/2, n, 3n/2) as a function of the signal SPP-setting a s . These fringes 
are shown together with the expected coincidence fringes as predicted by Eq. (0). The data 
show that the coincidence-count rate as a function of the orientation of the signal SPP 
has a parabola-like shape. Additionally, when the idler SPP is re-oriented, we see that 
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the fringe shifts accordingly. Consequently, we find, according to the quantum-mechanical 
prediction, that the parabolic coincidence fringe depends only on the relative orientation of 
the SPPs. Analytical calculations in which one assumes a completely separable (as opposed 




FIG. 4: The measured coincidence fringes (circles and broken line segments to guide the eye) and 
the theoretical fringes (curves) as a function of the orientation a s of the SPP in the signal arm. 
The parabolic fringe shifts concordantly with the orientation on of the SPP in the idler arm. 

to entangled) two-photon state while OAM is conserved per pair, result in fringes of zero 
visibility. In the context of a quantum-mechanical description of the SPDC process, the 
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experimental data show thus entanglement. However, based on the present set of data, we 
cannot prove quantum non-locality using the two-dimensional Clauser-Horne-Shimony-Holt 
version of the Bell inequality , in contrast to what we predicted earlier (f| . 

The dimensionality of the Hilbert space in which each photon of the pair lives, must 
be D > 2 as can be concluded from the parabolic shape of the fringes in Fig. after all, 
for D = 2, as in polarization entanglement, a sinusoidal fringe would appear. Regretfully, 
the fringes do not directly betray the precise dimensionality; nevertheless, we can make an 
estimate based on the following considerations. 

The dimensionality D of the subspace probed in the experiment is either limited by the 



number o 
K 
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entangled modes emitted by the BBO crystal, expressed by the Schmidt number 



M, or by the (spatial) tnoda! bandwidth of the ^alyzers, A, To estimate K, we 




follow Law and Eberly's approximate analytical theory |l6j |. 

Our experimental values for the pump-beam waist wq, crystal length L and the pump 
wavenumber k p result in a lower limit to the Schmidt number, K > 3700 ± 100. A com- 
parison between the approximate and exact theories for the value of K shows that, for high 
values of K, as in the present case, the approximate theory yields a severe underestimation. 

To make sure that our analyzers are not limiting the dimensionality of the probed sub- 
space to D < 3700, we need to establish that their modal bandwidth N ^> 3700. This can 
be accomplished by examining the conversion efficiency of the two analyzers, positioned as 
an optical train in a classical-optics experiment: we use one analyzer to convert the funda- 
mental Gaussian mode from a laser beam to an extended superposition of LG modes, while 
the other is used to convert this superposition back to the fundamental Gaussian mode. 
The conversion efficiency, which we found to be 98%, basically tells us how well the modal 
spectra of the two analyzers overlap. Combined with knowledge of the modal spectrum of a 
single analyzer (cf. Fig. |2J), this value yields the modal bandwidth iV of the analyzers. The 
experimental value of 98% [9], which is the fidelity of one analyzer mode as measured by 
the other analyzer, corresponds to a modal bandwidth N « 125000 (which corresponds, in 
fact, to the diffraction limit of the full optical setup). Consequently, the analyzers are not a 
bottleneck in this experiment and we can therefore conclude that the dimensionality of the 
probed Hilbert space is 3700 < D < 125000. 



In conclusion, we have reported an experiment on SPDC twin-photons demonstrating, 
within the quantum-mechanical description of SPDC, the entanglement of twin photons 
with respect to their fractional OAM. As this experiment involves a unitary transformation 
of basis from the natural LG states with integer OAM to a new basis with fractional OAM, 
the data explicitly shows that a very large part of the spatial Hilbert space is entangled, 
the dimensionality of which we estimate as D > 3700 per photon. The entangled space 
that we have probed in this experiment only represents one particular cross section of the 
full Hilbert space of (spatial) two-photon states. Presently, we are designing analyzers with 
alternative spatial geometries and step heights, and, additionally, analyzers that operate 
on both amplitude and phase of the photon state. Experiments with such analyzers will 
possibly allow us to sample a large variety of cross sections of the full Hilbert space. 
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